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In  this  paper,  from  the  viewpoints  of  both  the  first  and  the  second  law  of  thermodynamics, 
we  conduct  a  comprehensive  study  on  hydrogen-air  counter-flow  diffusion  combustion  in 
various  modes.  The  effects  of  air  inlet  temperature  (Tox;)  and  effective  equivalence  ratio  of 
fuel  (cp)  on  the  reaction  zone  structure  and  entropy  generation  of  combustion  are  revealed 
over  a  wide  range  of  Toxj  and  cp.  Through  the  present  work,  five  interesting  features  of 
combustion  of  this  kind,  which  are  quite  different  from  that  reported  in  the  literature,  are 
presented.  Especially,  for  the  first  time  we  divide  various  combustion  modes  in  the  cp  -  Tox; 
map  instead  of  the  popular  way  used  in  previous  studies.  Such  innovation  can  help  judge 
the  final  combustion  regime  more  straightforwardly  for  any  given  operative  condition. 
Copyright  ©  2011,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights 

reserved. 


1.  Introduction 

Hydrogen-air  counter-flow  diffusion  combustion  has  been 
extensively  studied  experimentally  or  computationally  in  the 
combustion  community  for  almost  40  years  since  it  provides 
a  nontrivial  research  prototype  to  deepen  our  insight  in 
detailed  chemical  kinetics,  pollutant  production,  ignition- 
extinction  response  of  flames  governed  by  high-temperature 
chemistry,  turbulence-combustion  interaction  and  so  on 
[1—4].  Especially,  the  recent  dilemma  caused  by  the  sharp 
increasing  demand  on  electricity  and  the  urgency  to  reduce 


man-made  emissions  of  green  house  gases  re-highlights  the 
importance  of  investigating  hydrogen-air  combustion  since 
hydrogen  is  one  of  the  most  promising  alternatives  to 
conventional  fossil  fuels  in  the  near  future  [5—10],  In  order  to 
bridge  the  gap  between  the  apparently  contrasting  goals  of 
minimization  of  pollutant  emissions  and  improvement  of 
process  efficiency,  some  novel  combustion  technologies,  such 
as  the  so-called  “flameless  oxidation”  [11]  or  “mild  combus¬ 
tion”  [12,13],  have  been  adopted  for  hydrogen-air  reactive 
flows  [14-19],  In  Ref.  [14],  Park  and  his  cooperators  numeri¬ 
cally  studied  the  flame  structure  in  mild  combustion  regimes 
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Nomenclature 

y 

relative  entropy  generation  rate 

A 

aspect  ratio 

Subscripts  and  superscripts 

a 

strain  rate 

0 

reference 

D 

diffusivity 

total 

total 

F/A 

hydrogen-air  ratio 

air 

air 

L 

half  distance  between  jets 

max 

maiximum 

P 

pressure 

si 

self-ignition 

Re 

Reynolds  number 

fuel 

fuel 

S 

local  entropy  generation  number 

oxi 

oxidizer 

T 

temperature 

in 

inlet 

u 

fluid  velocity 

St 

stoichiometric  value 

w 

width  of  jet 

c 

critical  value 

Greek 

AT 

V 

<p 

symbols 

temperature  difference 
kinematic  viscosity 
effective  equivalence  ratio 

cond 

chem 

mix 

vis 

heat  transfer 

chemical  reaction 

mixing 

fluid  friction 

of  hydrogen-air  laminar  flames  diluted  with  steam  and  found 
that  there  existed  an  oxidizer-side  temperature  limit  that  the 
combustion  mode  would  change  from  high  temperature 
combustion  to  mild  combustion.  Later,  Mollica  et  al.  [15] 
revealed  the  effect  of  preheating,  further  dilution  provided 
by  inner  recirculation  and  of  radiation  model  for  a  laboratory- 
scale  hydrogen-air  mild  burner  with  the  aid  of  CFD  technol¬ 
ogies.  In  Ref.  [16],  the  authors  found  that  pure  hydrogen  could 
not  reduce  thermal  NOx  emission  in  the  flameless  combustion 
regime.  With  the  aid  of  numerical  experiments,  Mardani  et  al. 
[17]  demonstrated  that  molecular  diffusion  could  not  be 
neglected  under  the  mild  combustion  condition  and  the 
influence  of  molecular  diffusion  on  mild  combustion  would 
increase  with  increasing  H2  in  the  fuel  blends.  Recently, 
Mardani  and  Tabejamaat  [18]  and  Wang  et  al.  [19]  investigated 
the  effects  of  hydrogen  addition  on  the  reaction  zone  struc¬ 
ture  of  a  JHC  burner  in  the  mild  regime.  Their  results  both 
showed  that  hydrogen  addition  could  organize  mild 
combustion  more  easily.  However,  as  clearly  shown  in  the 
above  literature  survey,  the  existing  studies  almost  all  are 
limited  in  co-flow  combustion  and  a  necessary  comparison  of 
hydrogen-air  counter-flow  diffusion  combustions  in  various 
combustion  modes  is  absent  yet.  Moreover,  they  are  based  on 
the  first-law  of  thermodynamics  and  ignored  the  more 
complicated  analysis  for  entropy  generation,  although  it  has 
been  widely  accepted  that  the  latter  is  better  for  optimum 
design  [20-23].  The  pioneering  studies  of  the  second  law 
analysis  on  opposing  jets  counter-flow  combustion  were 
conducted  by  the  present  authors  [6,24—26],  In  Refs.  [6,24—26] 
it  was  observed  that  the  characteristics  of  entropy  production 
in  opposing  jets  counter-flow  combustion  are  quite  different 
from  that  in  co-flow-jet  combustion.  Howerver,  the  discus¬ 
sions  in  Refs.  [6,24-26]  are  mainly  limited  in  conventional  air 
combustion  mode  and  a  systematical  investigation  on 
hydrogen-air  counter-flow  diffusion  combustions  in  various 
regimes  based  on  the  second  law  analysis  still  remains  as  an 
unsettled  question. 

The  main  originality  of  the  present  work  is  to  systemati¬ 
cally  investigate  hydrogen-air  counter-flow  diffusion 
combustion  in  various  combustion  modes  from  the 


viewpoints  of  both  the  first  and  second  law  analyses.  The  goal 
of  this  paper  is  to  reveal  the  influences  of  air  inlet  temperature 
(T0xi)  and  effective  equivalence  ratio  of  reactants  (<p)  on  the 
reaction  zone  structure  and  irreversibilities  generation.  The 
above  literature  survey  shows  no  open  study  on  this  topic. 
Especially,  in  this  paper  the  different  combustion  modes  are 
presented  in  a  q>  -  T0Xi  map  for  the  first  time,  which  provides 
a  more  practical  and  easier  way  than  that  in  previous  studies 
[12,14,31]  for  combustion  specialists  and  engineers  to  identify 
final  combustion  regimes.  The  governing  equations  for  flow 
and  scalar  fields  as  well  as  part  of  the  entropy  generation 
equation  are  solved  by  the  lattice  Boltzmann  (LB)  method 
recently  developed  by  the  present  authors  [27—29]  instead  of 
traditional  numerical  methods.  As  shown  in  our  previous 
studies  [26,29],  the  LB  method  possesses  significant  advan¬ 
tages  over  traditional  CFD  methods  for  the  analyses:  first,  the 
LB  method  is  more  suitable  for  massive  parallel  computing  of 
engineering  flows  and  more  efficient  to  evaluate  the  entropy 
generation  due  to  fluid  friction;  second,  the  LB  method  can 
accurately  represent  microscopic  physics  such  as  molecular 
diffusion  due  to  its  kinetic  nature,  which  can  be  neglected  in 
a  conventional  combustion  regime  [30]. 


2.  Specification  of  the  problem  and 
mathematical  modeling 

In  the  early  studies,  generally  the  mild  combustion  is  imple¬ 
mented  by  the  means  of  Hot-Oxidant-Diluted-Oxidant  [12]. 
Later,  Cavaliere  et  al.  organized  the  mild  combustion 
successfully  by  the  Hot-Oxidant-Diluted-Fuel  [31]  and  Hot- 
Fuel-Diluted-Fuel  [32]  fed  configurations.  According  to  the 
original  classification  by  Cavaliere  and  de  Joannon  [12]  under 
the  Hot-Oxidant-Diluted-Oxidant  fed  condition,  the  combus¬ 
tion  modes  can  fall  into  three  categories  except  “no  combus¬ 
tion”:  feedback  combustion,  high  temperature  air  combustion 
(HiTAC)  and  mild  combustion  according  to  the  difference 
between  the  inlet  temperature  Tin  (Tjn  =  TOXjfor  the  Hot- 
Oxidant-Diluted-Oxidant  and  Hot-Oxidant-Diluted-Fuel 
conditions  while  Tin  =  Tfuel  for  the  Hot-Fuel-Diluted-Fuel 
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configuration)  and  the  temperature  rise  AT  (AT  =  Tmax  -  Tin, 
where  Tmax  is  the  maximum  temperature  occurring  in  the 
reactor  and  Tsi  is  the  self-ignition  temperature  of  reactant 
mixture),  as  shown  in  Table  1.  The  modes  under  the  Hot- 
Oxidant-Diluted-Fuel  condition  are  similar  to  those  of  Hot- 
Oxidant-Diluted-Oxidant  configuration  [31].  However,  in  the 
Hot-Fuel-Diluted-Fuel  fed  operation,  Cavaliere  et  al.  [32]  found 
that  between  the  feed-back  combustion  region  and  the  no¬ 
combustion  region  there  would  be  an  additionally  transi¬ 
tional  zone  where  a  flame  thickening  occurred  due  to  heat 
release  but  pyrolytic  regions  disappeared.  This  transitional 
zone,  named  “flameless”  by  Cavaliere  et  al.  [32],  can  be 
roughly  characterized  by  Tfuei  <  Tsi  and  AT  >  0. 

The  configuration  of  planar  opposing  jets  illustrated  in 
previous  studies  [31,33],  namely  the  Hot-Oxidant-Diluted-Fuel 
fed  fashion,  is  adopted  in  this  study  because:  firstly,  various 
combustion  regions  all  can  be  implemented  by  it  with 
a  straightforward  way  [6,24-26];  secondly,  the  oxygen 
concentration  in  air  flow  is  kept  as  a  constant  (21%  by  volume 
ratio  in  the  present  study)  under  the  Hot-Oxidant-Diluted-Fuel 
fed  condition,  so  the  corresponding  analysis  can  be  simplified. 
It  is  well-known  that  the  influence  of  variable  oxygen  concen¬ 
tration  on  reaction  zone  structures  is  extremely  complicated 
[30]  and  accordingly  it  makes  a  comparison  in  a  consistent 
manner  very  difficultly.  For  example,  the  self-ignition  temper¬ 
ature  of  reactant  mixture  in  the  Hot-Oxidant-Diluted-Oxidant 
fed  mode  will  change  significantly  with  diluted  oxidizer 
concentration  and  HiTAC  cannot  be  investigated  under  the 
Hot-Fuel-Diluted-Fuel  fed  condition.  However,  these  deficits 
can  be  avoided  by  the  Hot-Oxidant-Diluted-Fuel  fed  fashion. 

The  problem  domain  and  boundary  conditions  are 
summarized  in  Fig.  1.  Two-dimensional  rectangular  coordi¬ 
nates  are  used.  The  origin  of  coordinates  is  located  at  the 
domain  geometric  center.  Two  parallel  stationary  walls  are 
located  at  y  =  ±L.  The  aspect  ratio  A  =  (L/W )  =  0.6.  The 
hydrogen,  diluted  by  nitrogen,  is  uniformly  ejected  from  the 
bottom  wall  with  temperature  Tjuei  =  T0  and  the  air,  consisting 
of  oxygen  and  nitrogen,  is  uniformly  ejected  from  the  top  wall 
with  temperature  Tair.  The  counter  flow  impacts  and  reacts  in 
the  reaction  zone.  Then,  the  diffusion  stagnation  “flame”  is 
formed.  The  burned  gas  flows  outward  along  the  x  -direction. 

There  is  one  important  parameter  for  the  combustible 
mixture:  the  effective  equivalence  ratio  </>,  which  is  given  by 
[5,24] 


F/A 

*~W. M 


(i) 


where  F/A  is  the  hydrogen-air  ratio  and  (F/A)st  refers  to  the 
stoichiometric  value  of  F/A  [24]. 

The  thermodynamic  and  transport  properties  appearing  in 
the  governing  equations  are  given  in  Ref.  [34].  The  detailed 


Table  1  -  The  classification  of  different  combustion 
modes  in  [12]. 

Combustion  mode 

Inlet  conditions 

Working  conditions 

Feedback  combustion 

Tin  <  Tsi 

AT  >  Tsi 

HiTAC 

Tin  >  Tsj 

AT  >  Tsi 

Mild  combustion 

Tin  >  Tsi 

AT  <  Tsi 

Fig.  1  -  Schematic  configuration  and  coordinate  system  of 
the  computational  domain. 


hydrogen-air  chemical  reaction  mechanisms  also  can  be 
found  in  Ref.  [34],  In  the  present  study,  uair  =  0.1  is  the  mean 
dimensionless  air  inlet  velocity;  p0  =  1/3  and  T0  =  1  are  the 
dimensionless  environmental  pressure  (1  bar)  and  tempera¬ 
ture  (300  K)  respectively.  The  dimensionless  characteristic 
length  is  L  =  1.  The  inlet  Reynolds  number  is  defined  as  [28] 

Re  =  U=L,  (2) 

vair 

where  vair  is  the  kinematic  viscosity  of  air.  The  velocity  ujuei  is 
determined  by 

Re  =  ^  (3) 

Vfue  I 

where  j^ei  is  the  kinematic  viscosity  of  fuel  mixture. 

The  mathematical  model  is  presented  in  two  parts.  The 
first  part  deals  with  the  simulation  of  stagnation  flames 
confined  by  planar  opposing  jets,  based  on  the  solution  of  the 
conservation  equations  for  reacting  counter-flow.  The 
dimensionless  governing  equations  in  Cartesian  coordinates 
for  such  laminar  steady  reacting  flows  have  been  given  in  our 
previous  studies  [27,28]  and  a  simple  LB  model  was  proposed 
in  them  to  solve  the  governing  equations.  The  second  part 
deals  with  the  computation  of  irreversibilities  or  entropy 
generation  in  counter-flow  combustion,  and  the  detail  on  it 
can  be  found  in  our  previous  publications  [6,24-26],  What 
should  be  mentioned  is  that  in  the  present  study  all  quantities 
used  in  computation  are  dimensionless.  The  detailed 
normalized  process  can  be  found  in  Refs.  [28,29]. 


3.  Results  and  discussions 

In  the  present  study,  we  investigate  the  effects  of  the  effective 
equivalence  ratio  <p  and  air  inlet  temperature  Toxi,  which  are 
important  operative  parameters  for  counter-flow  combustion, 
on  the  reaction  zones  and  irreversibilities  generation.  The 
variable  range  of  <p  is  very  wide,  from  the  fuel-ultra-lean 
region  (cp  <  0.7)  to  the  fuel-rich  region  (<p  =  2.0)  and  Toxi 
varies  from  T0  to  6  T0.  If  no  special  statement  made,  the  inlet 
Reynolds  number  is  fixed  by  Re  =  100.  Such  operating 
parameters  are  commonly  found  in  practical  applications  as 
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Fig.  2  -  Distributions  of  temperature  with  different  <p  and  Toxi  along  line  x  =  0:  (a)  Toxi  =  1  (b)  Toxi  =  3  (c)  Toxi  =  5  (d)  Toxi  =  6. 


well  as  fundamental  research  [1,24],  The  grid  resolution  used 
in  this  work  is  300  x  180.  It  has  been  demonstrated  in  our 
previous  works  [6,24,25,27,28]  that  such  grid  resolution  is  fine 
enough  to  obtain  grid-independent  numerical  solutions  for 
the  reactive  flow  in  this  study.  The  numerical  code  used  here 
is  described  and  validated  in  more  detail  in  Refs.  [6,24-26], 

3.1.  Reaction  zone  structures  in  different  combustion 
modes 

Following  the  way  in  Refs.  [31,33],  the  temperature  distribu¬ 
tions  are  illustrated  in  Fig.  2  to  reflect  the  characteristics  of  the 


reaction  zone  structure  within  different  combustion  modes, 
in  which  <p  =  0.5  (ultra-lean/higly  diluted  fuel  region), 
ip  =  0.7(critical  point  of  ultra-lean/lean  fuel  region),  <p  =  1 
(stoichiometric  reaction)  and  <p  =  1.2,  2  (rich  fuel  region)  with 
Toxi  =  1,3,5  and  6  are  chosen  as  the  representatives  for  the 
cases  investigated  in  the  present  work.  The  self-ignition 
temperature  of  hydrogen  in  air  at  the  present  strain  rate  is 
about  858  K,  corresponding  to  Tsi  =  2.86,  so  for  Toxj  <  Tsi  a  high 
temperature  source  is  used  to  ignite  the  reactants  [28].  As 
demonstrated  in  Fig.  2,  it  is  clear  that  without  preheating  the 
air  flow,  the  combustion  cannot  take  place  at  ip  <  0.5.  What¬ 
ever  Toxi  is,  the  temperature  peak  increases  with  tp  but  the 


Fig.  3  -  Variations  of  (a)  maximum  temperature  of  reactants  (b)  temperature  rise  with  inlet  air  temperature  according  to 
effective  equivalence  ratio. 
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Feedback  combustion 

HiTAC 

3.5 

Tsi  3.0- 

Mild  combustion 

2.5 

4- 

/  HiTAC 

<  2.0- 

"Flameless" 

3- 

1.5- 

Mild  combustion 

1.0- 
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NO  combustion  \ 
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<P 


Fig.  4  -  Combustion  modes  presented  in  (a)  TOXi  —  AT  map  (b)  <p  —  TOJ[i  map  at  Re  =  100. 


increment  of  flame  temperature  decreases  against  Toxi.  Since 
T0Xi  >  5,  the  temperature  peak  of  ip  =  0.5  becomes  vanishing, 
showing  the  special  feature  of  flameless  combustion 
mentioned  in  Ref.  [33].  These  observations  are  illustrated 
more  clearly  in  Fig.  3.  In  Ref.  [14],  it  was  observed  that  the 
thickness  of  reaction  zones  of  steam-diluted  hydrogen-air 
counter-flow  diffusion  combustion  increases  with  Toxi.  Their 
conclusion  can  be  confirmed  by  the  present  work,  although 
nitrogen  instead  of  steam  is  used  as  dilution  in  this  study. 
According  to  the  same  figure,  it  is  also  obvious  that  the 
thickness  of  reaction  zones  increases  with  <p ,  which  results 
from  that  the  reaction  zone  thickness  is  related  to  v/( D/a ) 
where  D  is  a  diffusivity  and  a  a  strain  rate  [14],  For  a  given 
strain  rate  (determined  by  the  inlet  Reynolds  number  [6]),  the 
mixture  diffusivity  increases  with  <p  since  the  diffusivity  of 
hydrogen  is  higher  than  nitrogen,  etc. 

Fig.  3  plots  the  maximum  temperature  of  reactants  Tmax 
and  temperature  rise  AT.  As  shown  in  this  figure,  Tmax  and  AT 
both  are  monotonic  linear  increasing  functions  of  <p.  Since 
ip  >  0.6,  these  lines  are  almost  parallel  with  each  other.  Similar 
observations  have  been  reported  in  Ref.  [14]  but  their  discus¬ 
sions  were  limited  in  mild  and  HiTAC  modes  and  hydrogen 
was  diluted  by  steam  instead  of  nitrogen.  Through  the  present 
study,  it  can  be  confirmed  that  these  characteristics  go 
throughout  the  entire  combustion  regions  of  hydrogen-air 
counter-flow  diffusion  “flame”  and  no  matter  which  dilution 


Fig.  5  -  Distribution  of  temperature  with  cp  =  0.4  and 
T0xi  =  2  along  line  x  =  0. 


(nitrogen  or  steam)  is  used.  As  cp  <  0.6,  because  the  heat 
release  due  to  chemical  reaction  decreases  significantly  with  <p 
and  becomes  so  slight  that  no  temperature  peak  appearing 
(namely  AT  ->  0),  therefore  the  curves  become  overlapping 
with  each  other. 

According  to  Fig.  3,  the  map  spanned  by  Toxi  -  AT  is  given  in 
Fig.  4  (a).  In  existing  open  publications  [12,14,31],  different 
combustion  modes  are  identified  with  the  aid  of  this  type  of 
map.  In  Ref.  [31],  Cavaliere  et  al.  identified  the  borders 
between  different  combustion  modes  for  methane-air 
counter-flow  diffusion  flame  under  the  Hot-Oxidant-Diluted- 
Fuel  condition.  They  found  that  the  cases  investigated  by 
them  fell  into  four  regimes:  NO  combustion,  feedback 
combustion,  HiTAC  and  mild  combustion.  However,  for 
hydrogen-air  counter-flow  diffusion  flame  under  the  Hot- 
Oxidant-Diluted-Fuel  condition,  besides  the  four  combustion 
modes  mentioned  by  Cavaliere  et  al.,  we  observe  an  addi¬ 
tionally  transitional  region  between  the  NO  combustion 
regime  and  the  feedback  combustion  regime,  as  shown  in 
Fig.  4(a).  Within  this  transitional  region,  although  the  pre¬ 
heated  air  temperature  is  not  high  enough  to  ignite  the  reac¬ 
tants  or  cause  fast  combustion  reactions,  the  accumulation  of 
the  input  heat  energy  carried  by  preheated  air  and  the  heat 
release  by  exothermic  reaction  can  sustain  slow  or  mild 
combustion  reactions.  The  case  with  rp  —  0.4  and  Toxi  =  2.0  is 
taken  as  the  example  to  explain  it  more  clearly,  as  illustrated 
in  Fig.  5.  In  this  figure,  it  can  be  seen  that  the  peak  temperature 
of  reactants  is  slightly  higher  than  Tsi  so  the  diffusion 
combustion  can  reach  steady  instead  of  extinction.  As 
mentioned  above,  the  flame  will  be  extinct  for  <p  =  0.4  if 
T0Xi  =  1-0  owing  to  the  heat  release  under  so  highly  diluted  fuel 
condition  is  too  small  to  sustain  combustion.  This  transitional 


Table  2  -  The  classification  of  different  combustion 
modes  for  hydrogen-air  counter-flow  under  Hot-Oxidant- 
Diluted-Fuel  condition. 


Combustion  mode 

Inlet  conditions 

Working  conditions 

Feedback  combustion 

Toxi  <  Ts; 

AT  >  Tsi 

HiTAC 

Toxi  >  Tsi 

AT  >  Tsi 

Mild  combustion 

T  •  >  T  • 

1  OXl  ^  1  SI 

AT  <  Tsi 

“Flameless” 

T  ■  <  T  ■ 

1  OX I  ^  1  SI 

AT  >  Tsi  -  T„xi 
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Fig.  6  -  Distributions  of  mass  fraction  of  H2)  02  and  H20  with  (a)  <p  =  0.4  (b)  <p  =  0.7  (c)  ip  =  1  (d)  ip  =  2  at  TOXi  =  3. 


zone  can  be  characterized  by  Toxi  <  Tsj  and  AT  >  Tsi  -Toxi. 
Consequently,  for  hydrogen-air  counter-flow  diffusion 
combustion  under  the  Hot-Oxidant-Diluted-Fuel  condition, 
there  exist  five  modes,  as  summarized  in  Table  2  and  the 
Toxi  -  AT  map  for  them  is  quite  different  from  that  in  Ref.  [31]. 
In  fact,  Cavaliere  et  al.  [31]  have  assumed  that  the  classifica¬ 
tion  of  various  combustion  modes  might  rely  on  fuels  used. 
The  present  results  confirm  their  assumption:  the  possible 
combustion  mode  for  hydrogen  is  quite  different  from  that  for 
methane.  It  is  very  interesting  that  the  transitional  region 
observed  in  the  present  work  is  similar  with  that  reported  in 
Ref.  [32].  In  Ref.  [32],  the  authors  also  found  that  there  was 
a  transitional  region,  being  roughly  characterized  by  1/ueI  <  fsi 
and  AT  >  0,  between  NO  combustion  and  feedback  combus¬ 
tion  regimes  and  named  it  “flameless”.  Accordingly  we  also 
called  the  transitional  region  observed  in  the  present  study  as 
“flameless”  because  the  major  features  of  these  two  transi¬ 
tional  regions  are  similar.  Nevertheless,  bear  in  mind  that  in 
Ref.  [32]  the  Hot-Fuel-Diluted-Fuel  fed  fashion  was  adopted  for 
methane-air  counter-flow  diffusion  combustion  while  in  the 
present  work  the  Hot-Oxidant-Diluted-Fuel  fed  condition  is 
used  for  hydrogen-air  reacting  flow. 

Originally,  the  identification  for  different  combustion 
modes  is  based  on  the  reactants  in  a  well-stirred-reactor  [12]. 
Soon  after  it  has  been  aware  in  the  combustion  community 
that  although  such  identifying  way  is  rigorous  but  it  needs  to  be 
evaluated  on  a  tassonomic  ground  [31].  Until  now  there  have 
been  several  attempts  made  along  this  direction,  for  example 
to  use  a  counter-flow  reactor  to  reveal  the  characteristics  and 
correspondingly  operative  conditions  in  different  combustion 


regimes  [31],  However,  more  efforts  are  desired.  For  instance, 
up  to  date  all  aforementioned  related  studies  discussed 
combustion  modes  in  the  Toxi  -  AT  map  like  that  illustrated  in 
Fig.  4(a).  The  shortcoming  to  adopt  such  means  is  obvious:  one 
cannot  determine  the  combustion  regime  a  priori  because  the 
operative  conditions  are  given  by  <p,  Toxi,  etc  instead  of  AT. 
Consequently  the  profile  presented  in  a  Toxi  -  AT  map  cannot 
straightforwardly  provide  necessary  information  for  the 
research  on  a  practical  system.  In  order  to  bridge  this  gap,  we 
divide  different  combustion  regimes  in  a  <p  -  Toxi  map  besides 
a  Toxi  -  AT  map,  as  shown  in  Fig.  4(b).  In  fact  the  dilution  degree 
of  reactants  can  be  also  characterized  by  <p  which  is  more 
popularly  used  in  the  combustion  community.  Through 


Fig.  7  -  Reaction  rates  of  hydrogen  with  different  <pat 
Toxi  =  3  along  line  x  =  0. 
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Fig.  4(b),  one  can  judge  straightforwardly  the  final  combustion 
regime  for  certain  given  operative  conditions  (<p,  Toxi  and  Re  for 
counter-flow  combustion),  while  it  is  hard,  at  least  not 
convenient,  to  make  such  estimation  in  a  Toxi  -  AT  map. 

The  distributions  of  mass  fraction  of  major  species  with 
different  <p  are  presented  in  Fig.  6.  It  is  obvious  that  in  an  ultra- 
lean/highly-diluted  fuel  region  (e.g.  q>  =  0.4)  the  reaction  zone 
structures  are  very  similar  to  the  so-called  Lindn’s  premixed 
flame  regime  [33]  in  spite  of  in  the  present  study  the  reactants 
being  separated  and  inter-diffuse:  02  and  H2  concentration 
distributions  cross  over,  which  implies  that  the  fuel  and 
oxidizer  are  not  consumed  at  the  thin  “flame”  region.  Conse¬ 
quently,  the  thin  flamelet  assumption  popularly  used  in 
turbulent  combustion  simulations  stands  only  when  the 
concentration  of  fule  is  high  enough  (e.g.  j  >  0.7  in  the  present 
work).  That  is,  H2  and  02  are  consumed  at  the  thin  flame 
region  and  the  reaction  rate  can  be  regarded  as  sufficiently 
high.  Similar  phenomena  have  also  been  reported  for  ultra¬ 
lean  methane-air  counter-flow  diffusion  combustion  in 
previous  work  [33]  but  unfortunately  has  attracted  little 
attention.  In  fact  it  is  an  important  issue  in  the  combustion 
community  since  current  flamelet  models  are  popularly  used 
for  combustion  simulation.  However,  the  present  study, 
together  with  our  recent  work  [35],  demonstrates  that  under 
highly- diluted  reactants  conditions,  either  for  fuels  or  for 


oxidizers,  the  validity  of  flamelet  assumption  is  questionable 
and  more  efforts  are  desired  on  this  topic.  It  can  be  seen  more 
clearly  through  Fig.  7,  in  which  the  reaction  rates  of  hydrogen 
with  different  <p  are  illustrated.  With  the  aid  of  Fig.  7,  one  can 
observe  that  the  peak  value  of  reaction  rate  of  H2  increases 
sharply  with  <p  but  at  the  same  time  the  range  of  reaction  zone 
decreases  against  q>  quickly.  These  observations  imply  that  in 
the  mild  regime  the  fast  chemistry  assumption  adopted  in 
traditional  combustion  simulation  becomes  inappropriate 
and  the  scale  separation  between  chemical  reaction  and  fluid 
convection  will  be  difficult. 

3.2.  Entropy  analysis  for  different  combustion  modes 

Fig.  8  illustrates  the  distributions  of  entropy  generation  S  in 
different  combustion  regimes.  From  this  figure,  it  is  obvious 
that  in  various  combustion  modes  the  patterns  of  S  are 
similar:  within  the  zone  near  the  “flame”  (referred  to  as  Zone  I 
in  our  previous  studies  [6,24-26]),  the  distributions  of  S  almost 
are  straight  lines  parallel  to  the  abscissa  and  entropy  gener¬ 
ates  intensively;  in  the  rest  room  outside  Zone  I  (referred  to  as 
Zone  II  in  Refs.  [6,24—26]),  the  distributions  of  S  form 
concentric  toroid-coil-like  curves  with  respect  to  the 
geometric  centers  of  jets  except  near  the  corners  and  when 
approaching  to  the  geometric  centers  of  jets,  entropy 


C 


d 


Fig.  8  -  Distributions  of  entropy  generation  at  different  modes:  (a)  HiTAC  (q>  =  2  and  Toxi  =  6)  (b)  Feedback  combustion  (cp  =  2 
and  Toxi  =  1)  (c)  Mild  combustion  (<p  =  0.7  and  Toxi  =  6)  (d)  “Flameless”  (<p  =  0.7  and  TOXi  =  1). 
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a  b 


Fig.  9  -  Profiles  of  S,  S„is,  Scond,  Smix  and  Schcm  along  line  x  =  0  at  different  modes:  (a)  HiTAC  (</>  =  2  and  Toxi  =  6)  (b)  Feedback 
combustion  (q>  =  2  and  Toxi  =  1)  (c)  Mild  combustion  (<p  =  0.7  and  Toxi  =  6)  (d)  “Flameless”  (<p  =  0.7  and  Toxi  =  1). 


generation  becomes  slighter.  As  explained  in  our  previous 
studies  [6,24-26],  this  kind  of  patterns  of  S  results  from  that 
entropy  generation  in  Zone  I  is  dominated  by  irreversibility 
due  to  chemical  reaction,  heat  transfer  and  mass  transfer 
while  in  Zone  II  it  is  predominated  by  irreversibility  due  to 
fluid  friction,  which  can  be  observed  more  clearly  from  Fig.  9: 
within  the  reaction  zones,  the  main  contributors  for  irre¬ 
versibility  are  Scond  and  Scflem-  The  peak  value  of  S„is  is  the  order 
of  ten,  much  less  than  those  of  Scond  and  Schem 
(0(  104)  -  O(105)),  so  the  fluctuation  of  Suis  along  x  =  0  can 
hardly  be  observed  in  Fig.  9.  Similar  phenomena  appear  for 
Smix  whose  peak  value  is  of  0(1O3).  Through  the  present  work, 


it  can  be  confirmed  that  such  type  of  entropy  generation 
patterns  is  one  of  the  intrinsic  characteristic  of  counter-flow 
diffusion  combustion,  no  matter  in  which  combustion 
regimes. 

In  a  recent  study  Soroudi  and  Ghafourian  [36]  investigated 
entropy  generation  in  HiTAC  and  mild  regimes  of  natural  gas 
counter-flow  diffusion  combustion  and  they  concluded  that 
the  intensity  of  irreversibility  generation,  namely  the 
maximum  value  of  entropy  generation  number  Smax,  will 
decrease  as  the  dilution  become  more  intense.  However,  we 
find  it  is  not  always  true  for  hydrogen-air  counter-flow  diffu¬ 
sion  combustion.  From  Fig.  10  (a)  one  can  clearly  observe  that 


Fig.  10  -  Variation  of  (a)  Smax  (b)  Stotai  with  different  <p  and  Toxi  along  line  x  =  0. 
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only  when  Toxi  <  5,  Smax  will  vary  like  a  monotonic  increasing 
function  of  cp.  While  Toxi  =  6,  Smax  will  firstly  increase  with  cp 
and  then  slightly  decrease  against  cp  since  cp  >  1.  In  other 
words,  only  when  Toxi  <  5,  the  conclusion  that  Smax  will 
decrease  as  the  dilution  become  more  intense  can  hold  water 
for  the  cases  investigated  in  the  present  work.  So,  besides  <p, 
T0Xi  also  plays  an  important  role  to  determine  the  intensity  of 
entropy  generation.  Surprisingly,  however,  it  was  neglected  in 
previous  work  [36]  and  a  detailed  comparison  between 
different  fuels  is  warranted  in  future.  Apparently,  for  a  given 
Toxi  since  the  effective  equivalence  ratio  beyond  a  critical 
value  <pc  (e.g.  cpc~  1.2  for  T0Xi  =  3,  <pc ~  1.5  for  Toxi  =  5  and  cpc~  1.6 
for  T0xi  =  6),  Smax  with  higher  Toxi  is  smaller  than  that  with 
lower  Toxi  (referred  to  as  Zone  III  in  this  paper).  Referring  to 
Fig.  4,  it  is  very  interesting  that  Zone  III  fully  overlaps  with  that 
of  HiTAC.  So  this  phenomenon  may  result  from  that  in  HiTAC 
regime,  the  irreversibility  due  to  heat  transfer  will  exceed  that 
due  to  chemical  reaction  and  become  the  dominant  contrib¬ 
utor  to  entropy  generation  within  the  reaction  zones,  as 
illustrated  in  Fig.  11.  In  this  figure,  the  relative  total  entropy 
generation  rate  due  to  heat  transfer  (ycond  =  ScmdttotaI/St„tai), 
chemical  reaction  (ychem  =  Schemtota, /Stotal ) ,  fluid  friction 
(Tuis  —  Suis.totai/Stotai)  and  mixing  ( 7 mrx  —  Sm[x.totai/St0iai )  are 
plotted  versus  T0Xj.  As  shown  in  Fig.  9,  in  HiTAC  mode,  the 
profile  of  entropy  generation  is  characterized  by  two  peaks  of 
Scond  separated  by  one  peak  of  SChem  and  accordingly  there 
appear  two  obvious  peaks  of  S.  While  in  other  combustion 
regimes,  there  is  only  one  obvious  peak  of  S.  Consequently 


Smax  in  HiTAC  regime  may  not  increase  with  more  fed  fuel  due 
to  the  area  with  intensive  entropy  generation  expanding 
much  more  rapidly. 

Fig.  10  (b)  plots  Stota!  with  different  <p  and  Toxi.  It  is 
straightforward  that  the  total  quantity  of  irreversibility 
increase  with  both  <p  and  Toxi  and  it  is  interesting  that  these 
curves  are  nearly  parallel  with  each  other.  However,  as 
mentioned  above,  the  shares  of  various  contributors  to 
entropy  generation  changes  significantly  with  different  <p  and 
Toxi-  Fig-  11  clearly  suggests  that  for  higher  dilution  with 
smaller  T0Xi,  the  reaction  zones  are  predominant  by  irrevers¬ 
ibility  due  to  chemical  reaction  while  instead  by  that  due  to 
heat  transfer  under  lower  dilution  and  higher  Toxi  conditions 
and  the  cross  point  of  ych em  and  ycond  moves  towards  lower  Toxi 
as  cp  increases.  As  shown  in  Fig.  2,  it  results  from  that  the 
temperature  gradients  become  much  steeper  as  cp  and  T0Xi 
increase.  Through  Fig.  11,  one  also  can  observe  that  the  vari¬ 
ations  of  different  contributors  to  irreversibility  are  extremely 
similar:  yChem  and  ymix  both  are  linearly  decreasing  functions  of 
T0xi  while  ycond  and  yvis  both  are  linearly  increasing  functions 
of  T0xi-  Moreover,  7uis  occupies  the  most  share  of  irreversibility 
while  Ymjx  is  the  smallest  one,  which  have  also  been  observed 
in  our  previous  studies  [6,24-26],  Accordingly,  we  can  make 
the  conclusion  that  such  observation  is  also  one  of  the 
intrinsic  characteristics  of  counter-flow  reactive  flow, 
regardless  of  the  combustion  mode. 

The  variations  of  relative  total  entropy  generation  rate 
versus  cp  are  illustrated  in  Fig.  12.  It  is  demonstrated  that  the 


T  T 
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Fig.  11  -  Variations  of  relative  total  entropy  generation  rate  versus  TOXi  at  (a)  cp  =  0.7  (b)  cp  =  1  (c)  <p  =  2. 
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Fig.  12  -  Variations  of  relative  total  entropy  generation  rate  versus  <p  at  (a)  Toxi  =  1  (b)  Toxi  =  3  (c)  Toxi  =  5  (d)  Toxi  =  6. 


patterns  of  different  contributors  to  irreversibility  spanned  by 
<p  are  similar,  too.  However,  as  Tox;  >  3,  except  ym;x  which 
becomes  negligible  since  ip  >  1.0,  yCond.  Tchem  and  yuis  all  are  not 
monotonic  functions  of  q> :  ychem  will  increase  with  <p  firstly  and 
then  decrease  against  tp  while  ycond  and  yuis  both  change 
inversely.  The  maximum  of  yChem  reaches  around  <p~0A.  With 
the  aid  of  Fig.  12,  one  can  observe  that  in  HiTAC  mode,  the 
irreversibility  due  to  heat  transfer  is  the  predominant 
contributor  within  the  reaction  zones,  agreeing  with  the 
conclusion  made  above.  While  in  the  “flameless”  regime,  the 
reaction  zones  are  dominated  by  irreversibility  due  to  chem¬ 
ical  reaction,  as  shown  in  the  same  figure.  Furthermore,  the 
feedback  and  mild  combustion  regimes  are  occupied  by  SChem 
and  Sconi  together:  generally,  in  these  two  combustion  modes, 
entropy  generation  due  to  chemical  reaction  dominates  the 
reaction  zones  under  lower  Toxi  and  intensive  dilution  condi¬ 
tions  but  irreversibility  due  to  heat  transfer  becomes  the 
leading  contributor  within  the  reaction  zones  for  higher  Toxi 
and  <p. 


4.  Conclusion 

This  study  has  investigated  for  the  first  time  the  effects  of  air 
inlet  temperature  (Toxj)  and  effective  equivalence  ratio  of 
reactants  ( <p )  on  the  reaction  zones  and  entropy  generation  in 
hydrogen-air  counter-flow  diffusion  combustion.  Especially, 
our  investigation  focuses  on  comparing  the  structures  of 


reaction  zones  and  entropy  generation  in  different  combus¬ 
tion  regimes.  In  order  to  achieve  this  goal,  a  wide  range  of 
parameters  is  chosen:  Toxi  varies  from  1  to  6  and  <p  from  0.1  to 
2.  Through  the  present  work,  five  important  features  of  reac¬ 
tion  zones  and  entropy  generation  in  this  kind  of  combustion, 
which  are  quite  different  from  that  reported  previously,  are 
revealed: 

1.  According  to  Ref.  [31],  four  different  combustion  modes  for 
methane-air  counter-flow  diffusion  “flame”  occur  under 
Hot-Oxidant-Diluted-Fuel  condition.  However,  for 
hydrogen-air  counter-flow  diffusion  “flame”  under  the 
same  feeding  condition,  we  observe  an  additional  transi¬ 
tional  combustion  mode  named  “flameless”,  which  implies 
the  classification  of  various  combustion  modes  may  closely 
rely  on  fuels. 

2.  Under  ultra-lean/highly-diluted  fuel  conditions  the  reac¬ 
tion  zone  structures  are  very  similar  to  the  so-called  Lilian's 
premixed  flame  regime  [33]  despite  the  reactants  being 
separated  and  inter- diffusing. 

3.  In  a  recent  study  [36]  it  was  concluded  that  the  intensity  of 
entropy  generation  in  HiTAC  and  mild  regimes  of  natural 
gas  counter-flow  diffusion  combustion  will  decrease  as  the 
dilution  becomes  more  intense.  However,  we  find  it  is  not 
always  true  for  hydrogen-air  counter-flow  diffusion 
combustion  and  that  Tox;  also  determines  the  variable  trend 
of  the  intensity  of  entropy  generation  in  HiTAC  and  mild 
regimes.  In  other  words,  for  certain  Toxi,  the  intensity  of 
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entropy  generation  may  increase  as  the  dilution  becomes 
more  intense. 

4.  In  HiTAC  regime,  the  area  with  intensive  entropy  genera¬ 
tion  will  expand  very  quickly  with  more  fed  fuel,  which 
causes  the  maximum  value  of  entropy  generation  number 
to  decrease  inversely. 

5.  In  HiTAC  mode,  the  irreversibility  due  to  heat  transfer  is  the 
predominant  contributor  within  the  reaction  zones.  While 
in  the  “flameless”  regime,  the  reaction  zones  are  dominated 
by  the  irreversibility  due  to  chemical  reaction.  Further¬ 
more,  the  feedback  and  mild  combustion  regimes  are 
occupied  by  them  together. 

Besides  the  above  findings,  the  present  work  has  also 
revealed  the  division  of  various  combustion  modes  in  the 
<p  -  Toxi  map  instead  of  the  Toxj  -  AT  map  popularly  used  in 
previous  publications.  The  advantage  of  such  innovation  is 
obvious:  one  can  judge  straightforwardly  the  final  combustion 
regime  for  given  operative  conditions  a  prior  with  the  aid  of 
the  <p  -  Toxi  map.  However,  it  is  hard  to  do  that  with  the 
Toxj  -  AT  map.  Furthermore,  a  more  detailed  comparison 
between  different  fuels  is  desired  in  future  study  since  we  find 
that  the  division  for  different  combustion  modes  and  the 
variation  of  intensity  of  entropy  generation  both  depend 
closely  on  fuels  used. 
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